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Abstract. We present here how the gravothermal or Antonov's instability, which was 
originally formulated in the microcanonical ensemble, is modified in the presence of a 
cosmological constant and in the canonical ensemble. In contrast to the microcanonical 
ensemble, there is a minimum, and not maximum, radius for which metastable states exist. 
In addition this critical radius is decreasing, and not increasing, with increasing cosmological 
constant. The minimum temperature for which metastable states exist is decreasing with 
increasing cosmological constant, while above some positive value of the cosmological constant, 
there appears a second critical temperature. For lower temperatures than the second critical 
temperature value, metastable states reappear, indicating a typical reentrant phase transition. 
The two critical temperatures merge when the cosmological density equals one half the mean 
density of the system. 



1. Introduction 

Antonov in 1962 [Ij discovered an instability that later became known as Antonov's instability 
or gravothermal catastrophe [2]. It signified the beginning of the study of thermodynamics 
of self-gravitating systems O IH [5l [6] and the nowadays rapidly developing field of statistical 
mechanics with long-range interactions [71 18] . 

In the original formulation of Antonov's instability [H |2], the system consists of particles 
(stars) bound by a spherical shell with insulating and perfectly refiecting walls. The particles 
interact only with Newtonian gravity (no cosmological constant or any relativistic effects are 
present) and the number of particles and the energy of the system are constant, i.e. the system is 
studied in the microcanonical ensemble. It was found that the entropy of the system has no global 
maximum and that for ER/GM^ < -0.335 no equilibrium states exist. For ER/GM'^ > -0.335 
there exist local entropy maxima, i.e. metastable states, only for density contrast below the 
critical value {po/PR)cr = 709, where po, pR are the central and edge densities, respectively. 

In gravity the ensembles are not equivalent. In the microcanonical ensemble, stable equilibria 
with negative specific heat do exist. This negative specific heat region is replaced by a phase 
transition in the canonical ensemble [3]. The Antonov system in the canonical ensemble (with 
no cosmological constant), is studied in detail by Chavanis [1]. He found that the free energy 
does not have a global minimum and equilibria exist only for temperatures greater than the 
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critical dimensionless temperature {RT /GM)cr = 0.40 which corresponds to dimensionless 
inverse temperature {GM (3 / R)cr = 2.52. These equihbria are stable only for density contrast 
less than the critical value {po/PR)cr = 32.1 (see Ref. [9], as well). 

We studied the gravothermal instability with a cosmological constant A in the microcanonical 
ensemble in Ref. pLO]. It was found that for a fixed energy, the maximum radius (which equals 
R^^ = {-0.335 /E)GM'^ with no A) is increasing with increasing A, while for A > there appears 
a second critical radius which is decreasing and the two radii merge at some marginal A. This 
was identified as a reentrant behaviour (see Ref. |1H [T^ [THl [TJ] for details on reentrant phase 
transitions in statistical systems with long-range interactions). In the present work, the system 
is placed in a heat bath and the walls are non-insulating, i.e. we study the canonical ensemble. In 
contrast to the microcanonical one, we find that the critical radius Rcr, for a fixed temperature 
and mass, is decreasing with increasing cosmological constant. In addition the instability occurs 
for R < Rcr in the canonical ensemble, in contrast to the microcanonical one. This is evident in 
case A = 0, since then there exist no equilibria if GM(3/R > 2.52 =^ R < GM (3/2.52 for a fixed 
temperature and mass. A reentrant phase transition appears in the canonical ensemble, when 
one examines the system with respect to the temperature for a fixed radius. There appear two 
critical temperatures for A > 0. For the intermediate temperature values, no thermodynamic 
equilibria are possible. 



2. Analysis and results 

In the Newtonian limit, the Poisson equation is modified [151 ll6j in the presence of a cosmological 
constant A as: 

S7^(j) = AnGp - SttGpa (1) 
with p\ = Thus, the gravitational potential can be decomposed to two parts 

4> = (t^N + 4>A (2) 



where </>Ar is the Newtonian and (/>a the cosmological potential, that are given by 

^N = -Gf T$\d^r' , = -^pAr^ (3) 
J \r — r'\ 3 

For convenience we call the Newtonian limit of de Sitter (pA > 0) and Anti-de Sitter {p\ < 0) 
spaces [17j just dS and AdS, respectively. 

We consider a self-gravitating gas of N particles with unity mass inside a spherical non- 
insulating shell inside a heat bath and restrict only to spherical symmetric configurations. 
We work in the mean field approximation, for which the A^-body distribution function is 
replaced by the one body distribution function /(r, -u). The Boltzmann entropy is defined as 
S = —k J f\\ifd^rd?p. The Helmholtz free energy is equal to F = E — TS. Equivalently one 
can work with the Massieu function [H [9] J = —F/T that gives 

J = S-^E (4) 

The maximization of J with constant T and the maximization of S with constant E (and 
constant M in both cases), with respect to perturbations 5 p., are the same [3] to first order in 
5p , i.e. give the same equilibria, described by the Maxwell-Boltzmann distribution function 

3 

i-V p{r)e-'^P-' (5) 



where 

p{r) = poe-'^(^-<^(°« (6) 

This is proved in [3] without a cosmological constant, and holds in the presence of A as well, 
with the difference being that (p does not satisfy Poisson equation now, but equation ([TJ. The 
two ensembles have the same equilibria. However, the second variations of J and S are different. 
Therefore, what is different in the two ensembles is the turning point at which an instability 
sets in. This causes a great qualitative difference for the two ensembles. 

Introducing the dimensionless variables y = f3{(j) — </'(0)), x = ry/AnGpo/S and A = 2p/^/po, 
and using equation ([6]), equation ([1]) becomes 

x"^ dx 

called the Emden-A equation. Let us call z = R^J AttGpqP the value of x at R. In order 
to generate the series of equilibria needed to study the stability of the system, the Emden-A 
equation has to be solved with initial conditions y(0) = y'(0) = 0, keeping M constant and for 
various values of the parameters pA, /3, po- This is a rather complicated problem, since, unlike 
A = case, while solving for various z, mass is not automatically preserved, because of the mass 
scale Ma = p^^nR^ that A introduces. A suitable A value has to be chosen at each z. We define 
the dimensionless mass 

2Ma SttpaR^ 2pA 

where p is the mean density of matter. Calling z = R\/4ttGpqP the value of x at R, equation ^ 
can also be written as m = 3B/Xz'^ where B = GMf5/R is the dimensionless inverse temperature. 
It can be calculated by integrating the Emden-A equation, to get: 

B{z) = zy\z) + \\z^ (9) 

We developed an algorithm to solve equation ([7]) for various values of A, z keeping m fixed. 
From equation (l8|) it is clear that solving for various fixed m can be interpreted as solving for 
various pA and/or R for a fixed M. 

We find that in dS case there exist multiple series of equilibria for a given cosmological 
constant, while in AdS case there is just one series likewise A = case. Studying the series 
/3(^) for both AdS and dS, we determined the critical values of radii and temperature beyond 
which there are no equilibria. These critical values are plotted in Figures [1] and [21 with respect 
to PA) where the unshaded region is the region of instability in both figures. 

In Figure [1] the radius denoted Ra is the minimum radius for which thermodynamic equilibria 
exist in the canonical ensemble. Ra is the radius of the homogeneous solution which exists only 
in dS and is defined by equation p = 2pA = const, i.e. Rh = (3M/(87rpA))^. In the shaded 
region of Figure [T] there can always be found metastable states. For AdS case (pA < 0), these 
have always monotonically decreasing density (po > pr) and suffer a transition to unstable 
equilibria at a critical density contrast {po/pR)cr whose value depends on pA- For dS case and 
for R < Rh, the metastable states suffer a transition to unstable equilibria likewise AdS, but the 
density is not necessarily monotonic, neither does po > pr always hold. For each pA there exists 
a tower of solutions with qualitative (and quantitative of course) different density functions. In 
the region R > Rh, all metastable states have po < pr. The ones, whose density is a monotonic 
function of R do not suffer any transition to unstable equilibria, while all others series have a 
turning point. 

In Figure [2] we see that in AdS the critical temperature is decreasing for increasing 
cosmological constant, while in dS there appears a reentrant phase transition. For small positive 
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Figure 1. The critical radius w.r.t. pA for 

fixed M, (3. 
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Figure 2. The critical temperature w.r.t. p\ 
for fixed M, R. 



values of pA, compared to the mean density p, the system cannot stabilize for low temperatures. 
It needs higher temperature values to obtain a significant pressure gradient to balance gravity. As 
PA is increasing, the cosmological repelling force increases, enforcing the total outward pushing 
force and hence, enabling the system to stabilize at lower temperatures. Point A, denotes the 
marginal value p^*" = p/4 for which a dynamically static equilibrium is allowed. This limiting 
configuration corresponds to a singular solution at which all mass is concentrated at the edge. 
It can easily be calculated: GM/2R^ = ^pf'^R =^ p^" = 3M/16ttR^ = p/4. For pA > pT"", 
solutions with T = (static dynamical equilibria) do exist and as pA increases the mass is 
allowed to be closer to the center. All thermodynamic equilibria below the lower critical branch 
(low temperatures) have increasing density (po < pn). For fixed R, M and a fixed pA with 
p/4 < Pa < p/2, as we increase the temperature beginning from zero, we pass from a region 
of metastable states with po < pR to a region of no equilibria, at a first critical temperature 
Ti- We have to increase the temperature even more to reach a new region of metastable states 
at a second critical temperature T2. For T > T2 the metastable states have decreasing density 
Po > P-R if < Rh- This is a reentrant phase transition similar to the ones found in other 
statistical systems with long-range interactions 
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